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Abstract 



At high temperatures or densities, hadronic matter shows different forms of critical 
behaviour: colour deconfinement, chiral symmetry restoration, and diquark condensation. 
>• ! I first discuss the conceptual basis of these phenomena and then consider the description 

of colour deconfinement in terms of symmetry breaking, through colour screening and as 
CS| ' percolation transition. 

O 

o 
o 

^ ■ 1. States of Matter 

Qh! Hadronic matter is endowed with an inherent density limit. The usual hadrons have 

^ I an intrinsic size — 1 fm, so that a hadron needs a space of volume Vh — (47r/3) to 

exist. Therefore 

nc ^ ~ 0.25 fm"^ ^ 1.5 no (1) 



^ i is the limiting density for such a medium; here uq = 0.17 fm~^ denotes the density of 

normal nuclear matter. In turn, this also leads to a limiting temperature for hadronic 
matter, 

T, ~ (l/n,) ~ 0.2 GeV, (2) 

as pointed out by Pomeranchuk almost fifty years ago ^ . Considering an ideal gas of res- 
onances, whose composition was based on a classical partitioning problem 0, Hagedorn 
found a very similar limit, which he proposed as the ultimate temperature of strongly 
interacting matter [Q. Dual resonance dynamics led to an essentially equivalent com- 
position law 1^. Soon afterwards Cabibbo and Parisi noted that the 'limit' more likely 
corresponded to a critical point, signalling the transition to a new state of matter, the 
quark-gluon plasma 0. What remains is the realization that, on geometric, combinatoric 
or dynamic grounds, hadron thermodynamics defines its own limit. 



1 Invited talk at the 3rd Catania Relativistic Ion Studies ORIS 2000, Acicastello/Italy, May 22 - 26, 
2000; at the XL Cracow School of Theoretical Physics, Zakopane/Poland, June 3 -11, 2000; and at the 
Conference on Strong and Electroweak Matter SEWM 2000, Marseille/France, June 13 - 17, 2000. 
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This limit can be approached in two ways. The obvious is to compress cold nuclear 
matter beyond no; but in relativistic thermodynamics, 'heating' mesonic matter leads to 
particle production and thus also increases the density. As a result, strongly interacting 
matter has a T—fis phase diagram, where T denotes the temperature and the chemical 
potential specifying the overall baryon number density. In the T — fiB plane, there must 
thus be a limiting curve for hadronic matter, beyond which the density is too great to 
allow the existence of hadrons. 

Since QCD defines hadrons as bound states of quarks, the general phase structure of 
strongly interacting matter is quite evident: for densities below ric, it consists of colourless 
hadrons, i.e., colour singlet bound states of three quarks or of a quark-antiquark pair. 
Above He, deconfinement leads to a medium consisting of coloured constituents. There 
are three possible forms for such constituents: 

• coloured massive quark-gluon states: constituent quarks; 

• coloured massive quark-quark states: diquarks; 

• coloured massless quarks and gluons: a quark-gluon plasma. 
What actually happens in the different regions of the T—fiB diagram? 

Before addressing this question, we have to consider how deconfinement can take place. 
The confining potential between a static quark and antiquark separated by a distance r 
has the idealized form 

V{r) ~ err, (3) 

where a ~ 0.8 GeV/fm is the string tension. The quarks inside a hadron are therefore 
confined: the hadron cannot be broken up into its coloured constituents, since this would 
require an infinite amount of energy. 

In a dense medium, however, there is another way to dissociate bound states. The 
presence of many other charges leads to charge screening, which reduces the range of the 
forces between charges. A well-known example is Debye screening, which suppresses the 
long-range part of the Coulomb potential between electric charges, 

V{r) = ^ ^ ^ e-^^ (4) 

where td = /^"^ is the Debye radius, defining the range of the force remaining effective 
between charges in the medium. When it becomes smaller than the atomic binding radius, 
an insulator consisting of charge-neutral atoms turns into a conducting plasma of unbound 
electric charges [H. In QCD, the corresponding effect leads to colour screening. 



V{r) = ar ^ or 



-fir 



fir 



(5) 



where fi ^ now defines the colour screening radius^. Deconfinement thus is the insulator- 
conductor transition of QCD, with colourless bound states as constituents below and 

^The difference in the form of the screening functions in Eq. and (||) is due to the different forms 
of the unscreened potentials Q . 
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coloured constituents above the deconfinement point. But what is the nature of the 
conducting phase here? 

One way to study that is to consider the effective quark mass. The input quark masses 
in the QCD Lagrangian are (for u and d quarks) almost zero, rrig ~ 0. In the confined 
phase, hadrons behave as if they consist of constituent quarks of mass mg, with m„ ~ Stuq 
and TUp ~ 2mQ, for nucleons and (non-Goldstone) mesons, respectively. Hence here the 
quarks manage to 'dress' themselves with gluons to acquire a mass rriQ ^ 0.3 GeV. 

At sufficiently high temperatures, thermal motion will presumably 'shake off' the 
dressing, so that somewhere in the course of the hadron-quark matter transition there 
will be an effective quark mass shift mg —>■ rrig. For vanishing rriq, the QCD Lagrangian 
is chirally symmetric; hence this chiral symmetry must be spontaneously broken in the 
confined phase and restored in the hot QGP. One therefore often refers to the shift in 
effective quark mass mg — > as chiral symmetry restoration. Such a shift is more general, 
however, and can occur as well for ^ 0, as shown by the shift of the effective electron 
mass between insulator and conductor. 

That leads to the next question: is there also a colour superconductor? At low temper- 
atures, collective effects of an electrically conducting medium can overcome the Coulomb 
repulsion between like charges and lead to a binding of electrons into doubly charged 
Cooper pairs. These, being bosons, can condensate to form a superconductor. In QCD, 
the conditions for creating a superconductor are in fact much more favorable. An attrac- 
tive local potential couples two triplet quark states to a bosonic anti-triplet diquark state, 
so that in QCD there is a dynamical basis for colour superconductivity through diquark 
condensation 0. 

In very recent years, the low temperature, high baryon density part of the QCD phase 
diagram has received much renewed attention, resulting in the prediction of different 
superconducting phases and several transitions 0. Although of great theoretical interest, 
this region is for the time being accessible neither to lattice studies nor to experiment. I 
shall therefore concentrate here on the high temperature, low baryon density region and 
refer to for a discussion of colour superconductivity. 

Taking into account what was said, a first guess of the QCD phase diagram leads to 
a four-phase structure of the generic form shown Fig. |l|; as noted, the 'diquark phase' 
shown there may well consist of several different superconducting phases. In any case, 
however, lattice QCD tells us that Fig. [l| is wrong: at /i = 0, deconfinement and chiral 
symmetry restoration coincide, so that there is no constituent quark phase. One of the 
main points I want to address here is why this is so. 

A second guess is shown in Fig. |^, with hadronic, diquark and QGP phases. As far as 
we know, this one may well be correct; for /i^ = 0, it is, as we shall see from the results 
provided by finite temperature lattice QCD. 

In the classical study of critical behaviour, the analyticity of the partition function 
Z(T,iJ,,V) in the thermodynamic limit of infinite volume, V oo, governs the phase 
structure, with critical points defined through the divergence of derivatives of Z{T,fi, V). 
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Figure 1: Four-phase diagram of strongly interacting matter. 
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Figure 2: Three-phase diagram of strongly interacting matter. 
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The occurrence of such singularities can be attributed to different underlying physics 
mechanisms. In statistical physics, spontaneous symmetry breaking, charge screening 
and cluster percolation have received particular attention. In the following sections, I 
shall consider these three 'mechanisms' in the context of statistical QCD. 

2. Symmetry Breaking 

In contrast to the statistical mechanics of condensed atomic matter, the phase struc- 
ture of strongly interacting matter can be determined ab initio from QCD as input dy- 
namics, without intermediate models. Given the QCD Lagrange density 

jCqcd = -F.^F'^" - i^iz^d, - g^A^ + m,)^, (6) 

where the first term describes the pure gluon sector and the second the quark-gluon 
interactions, one defines the partition function at temperature T as 

Z{T,V) = J dAdijj dip exp{-S{C,T)}. (7) 



Here A denotes the gluon and ijj the quark fields; the QCD action is given by 

S(C;T,V)= I d^x f^^^ dr C(A(x,t),Mx,t)), (8) 
Jv Jo 

as integral over the volume V of the system and a slice of thickness 1/T in the imaginary 
time r = ixQ. From Z{T, V) one then obtains the usual thermodynamic functions; the 
derivative with respect to T leads to the energy density e(T), that with respect to V the 
pressure P{T), and so on. Both the dynamic input theory and its associated thermody- 
namics are thus completely specified; the problem lies in the evaluation, for which one 
has to resort to the numerical simulation [|1^ of the lattice formulation |T2| . 



The conventional deconfinement probe in finite temperature lattice QCD is the expec- 
tation value L[T) of the Polyakov loop |[T^, |T^. Through 



L(T) ~ lim e-^^''^/^, (9) 

r^oo 

it is related to the potential V{r) coupling a static quark-antiquark pair. In the confined 
phase, this potential diverges as V{r) for r — oo (see Eq. (^)), while in the deconfined 
phase it converges to a finite value. It is found that L{T) = in the temperature range 
T < Tc, defining the confinement region, and that L{T) > for T > Tc, specifying the 
deconfinement region. 

Strictly speaking, V{r) ~ ar diverges for r — oo only in a theory with infinitely heavy 
quarks. In real QCD with light dynamical quarks, the string breaks when it becomes 
energetically more favourable to produce a quark-antiquark pair, i.e., when V{r) ~ 2mQ, 
where mq is the mass of a 'dressed' constituent quark and 2mQ the mass of a (non- 
Goldstone) meson. The quark of this newly produced pair combines with the original 
antiquark, the antiquark with the original quark, thus making two strings out of one. In 
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full QCD, the Polykov loop therefore does not vanish in the confining region, but only 
becomes exponentially small, 



L{T) ~ e-2'"«/'^ for T < T^. (10) 

For rriQ ~ 0.3 GeV and Tc ~ 0.15 GeV, this makes L{Tc) ~ 0.02 instead of zero. Never- 
theless, L{T) is now no longer a real order parameter, and deconfinement therefore seems 
to be not really defined as a critical phenomenon in QCD with light dynamical quarks. 
We shall return to this problem several times and show how it might be solved. 

From statistical mechanics, it is known that phase transitions are generally associated 
to symmetries of the system. Thus the Hamiltonian of the simplest spin theory, the Ising 
model, is 

n = -jJ2s^sj-BY.Si, Si = ±iyt, (ii) 

i<j i 

where the first sum runs only over nearest neighours on the lattice, J denotes the exchange 
energy between spins and B an external magnetic field. For S = 0, ?^ is invariant under 
the global Z2 symmetry of flipping all spins, Si — —Si V i. The thermodynamic states of 
this system share this symmetry for T >Tc, where now is the Curie point; as a result, 
the expectation value of the spin, the magnetization m{T), vanishes in this 'disordered' 
temperature region. Below T^, however, the system becomes ordered, the spins choose to 
align either up or down, making m{T) 7^ 0. Since 'up' or 'down' are equally likely, the 
symmetry as such is preserved; the actual state of the system, however, spontaneously 
breaks it by choosing one or the other. The magnetization transition of the Ising model 
thus corresponds to the spontaneous breaking of the inherent global Z2 symmetry of the 
Ising Hamiltonian. 

This line of argument applies directly also to deconfinement in pure SU{N) gauge 



theories [13, 14]. The corresponding Lagrangian in lattice QCD 



£sU(N)(f/^fcf/fc^t/+f/+) (12) 

depends on the products of four SU (N) matrices Uij on the links of the smallest closed 
loops of the lattice. It remains invariant under a global 'flip' e C SU{N) of 
all matrices associated to a spatial hyperplane, with U^^r znU^^t V x; here = 
exp{r(27rz/A^)} with r = 1, 2, A^. 

The Polyakov loop, on the other hand, does not remain invariant under such global 
Zn transformations, with 

L ~ <ReTr U^^^U^^r > ^ zn L. (13) 

It is thus the analogue of the magnetization of the Ising model, in the sense that it tests if 
the state of the system shares or spontaneously breaks the symmetry of the Lagrangian. 
This feature becomes particularly transparent for SU{2) gauge theory, where Z2 = ±1, 
so that the transformation thus just means flipping the sign of the Polyakov loop, L — > 
—L. In the temperature region in which L{T) = 0, i.e., in the confinement region, the 
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states are ZAr-symmetric, while for deconfinement, with L{T) > 0, the Zn symmetry is 
spontaneously broken. In other words, deconfinement in pure SU{N) gauge theory can 
be defined as the spontaneous breaking of a global Zn symmetry of the corresponding 
Lagrangian. 

The similarity between spin and gauge systems goes in fact much further [jl5| . In Fig. 
^, we compare schematically the temperature behaviour of the Polyakov loop L{T) and 
the magnetization m{T), together with that of the corresponding susceptibilities Xl{T) 
and Xm{T)- The latter measure the fiuctuations of the order parameters at the transition 
point and thus diverge there. For SU{2) gauge theory as well as for the Ising model, the 
transition is continuous, and so near the functional behaviour in the two cases can be 
written as 

L(T) ~ (T - T,)^^ T > T,; ~ |T - T,|-^^ (14) 

and 

m(T) ~ (T, - Tf-, T < Te; x^T) ~ |T - T,!"^™, (15) 

where (3 and 7 denote the critical exponents for the two transitions. While SU{N) gauge 
theories in general have a more complex structure than spin theories, their critical be- 
haviour becomes in fact identical: they belong to the same universality class [^, which 
means that (3^ = and 7^ = 7^. The confinement /deconfinement transition in SU{N) 
gauge theories is thus structurally the same as the disorder /order transition in spin theo- 
ries; both are based on the spontaneous breaking of a global Z^- symmetry of the under- 
lying dynamics. 




Figure 3: Schematic temperature dependence of the magnetization m(T) and the Polyakov 
loop L(T), together with the corresponding susceptibilities. 

The introduction of dynamical quarks in full QCD explicitly breaks this Zn symmetry; 
it effectively adds a term to the SU{N) action which is proportional to L: 

Sqcd ~ 5'su(N) + i^{mg)L, (16) 

where n{mq) —>■ for ruq —>■ 00. Comparing Eqs. (|T^) and (|TTp, we see that dynamical 
quarks in a sense play the role of an external field B in spin theory. Just as B aligns the 
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spins and prevents m{T) from ever completely vanishing, so does rriq result in a Polyakov 
loop which is always non-zero. We had seen above that another way of arriving at this 
conclusion is through string breaking. Hence there must be some implicit relation between 
the effective external field and the constituent quark mass determining the string breaking 
point. The effect of dynamical quarks on the Polyakov loop is shown schematically in Fig. 
^. We note in particular that even for nig = 0, there is 'almost critical' behaviour, with 
a sharp variation at a temperature considerably below the deconfinement temperature of 
pure SU (3) gauge theory What are the reasons for this behaviour? 



L(T,mq) 




mq=OI 


mq=oo 




T 



Figure 4: Schematic temperature dependence of the Polyakov loop in pure SU (3) gauge 
theory (m^ = oo) and in full QCD with two massless quark flavours (m^ = 0). 

From the string breaking picture, we expect the effective external field acting on the 
Polyakov loop as generalized spin to be inversely proportional to the constituent quark 
mass, B ~ l/rriQ [^, [T^, [T^. Thus deconfinement should occur when rriQ — * 0. 

For rriq = 0, the QCD Lagrangian is chirally symmetric; however, the state of the 
system under given conditions need not share this symmetry. The chiral condensate, 
K{T) = {ipijj) ~ itlq, provides an order parameter to probe if and when the chiral 
symmetry of the Lagrangian is spontaneously broken. It is found that 

K{T) ^ implying mg ^ ^ T < T^, (17) 

and 

K{T) = implying mg = ^ T > T^. (18) 

where MeV is the chiral symmetry restoration temperature. The functional behaviour 
of K and the corresponding fluctuation susceptibility xk, 

K{T) ~ (T^ - TY'<, T < T^; xk{T) ~ |r - TJ"^-, (19) 

is illustrated in Fig. |^. For T < T^, K is large and hence the effective external field 
B ~ ^/rriQ ~ 1/K{T) is small, so that the Polyakov loop is almost disordered, implying 
confinement-like behaviour. At T = T^, B suddenly becomes large; it aligns the Polyakov 
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Figure 5: Schematic temperature dependence of the chiral condensate K{T) in full QCD, 
together with the corresponding susceptibility. 



loops, implying the onset of deconfinement. We thus find that chiral symmetry restoration 
induces colour deconfinement. 

These considerations have immediate consequences which can be tested in finite tem- 
perature lattice QCD. The magnetization m(T, B) for non-vanishing external field B 
becomes an analytic function of T and B. In full QCD, we therefore assume the Polyakov 
loop for mg 7^ to be an analytic function of T and K. This leads to 



From this we obtain 



and 

for the Polyakov loop susceptibilities Xm Xt- Since the chiral susceptibilities Xm = 
{dK/dmq)T and Xt = idK/dT)^^ diverge at T = T^, relations {^/^ imply that the 
Polyakov loop susceptibilities must share this singular behaviour, with the same critical 
exponents. Present lattice calculations for full QCD are not yet precise enough to allow a 
conclusive determination of critical exponents. In Figs. ^ and ^ it is seen, however, that 
the increase of the chiral susceptibilities Xm Xt — is indeed accompanied 

by a similar increase in the Polyakov loop susceptibilities Xm Xt- 

In QCD with massless dynamical quarks, the chiral condensate K and the Polyakov 
loop L thus are analytically related; there is only one critical point T = T^, at which both 
K{T) and L{T) exhibit non-analytic behaviour. At T^, the quarks loose their effective 
mass and at the same time become unbound. 



dK. 



(20) 



Am 



dT 



dL 



dL\ f dK 



^ dnig ) rp \ dK j j,\ dm. 



1 / T 



dL 



'dK 



dK)^\dT 



(21) 
(22) 
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Figure 6: The chiral susceptibilities Xk Xk functions of the temperature variable 

K = 6/(7^. 
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Figure 7: The Polyakov loop susceptibilities and Xk as functions of the temperature 
variable k = Q/g"^. 



3. Colour Screening 

In pure SU{N) gauge theory, the potential between a static quark and antiquark 
increases linearly and unbounded with their separation r, V{r) ~ err. In full QCD, the 
string breaks when V{r) ~ mg, even in vacuum, i.e., at T = 0. We can interpret this 
by attributing a screening behaviour to the sea of virtual massless quark-antiquark pairs. 
Starting from Eq. (|lD, we thus have 



V{r, /i) = ar 



1 



-/ir 



a 



-fir 



fir 



(23) 



10 



where the last term includes Coulomb and transverse string effects. In Eq. (pSQ, the string 
tension a and the Coulomb coupling a are taken to be constants, with /i(T) temperature- 
dependent. First, we now want to determine the vacuum screening mass [i(T = 0). 

In the spectroscopy of heavy quarkonia, such as the J/ip or the T, the masses and 
widths of all bound states are determined by solving the Schrodinger equation 

n^i = Mi^i, (24) 

with the Hamiltonian 

n = 2m^b + V{r,fi = 0) (25) 

rrich 



given in terms of the potential (|23|). Here i specifies the cc or bb bound state under 
consideration, Mj its mass. By comparison for the results to quarkonium data, the four 
constants in the potential are determined, giving a = 0.192 GeV^, a = 0.471 and = 
1.32 GeV, rrLh = 4.75 GeV for the bare charm and bottom quark masses, respectively. 
The vacuum screening mass can now be obtained by comparing the gap between the open 
charm or bottom threshold, 2Md or 2Mb, and a given quarkonium state Mj, to its gap 
with respect to the infinite range potential, 

E^i,, = 2Md,b -M, = 2me,6 + - - M,. (26) 



The result 



= TTTTT-^ 7 (27) 

gives /i(T = 0) ~ 0.18 GeV for both charmonium and bottonium states. The fact that 
the large difference between rric and rrib plays no role here is an indication that the states 
are indeed heavy enough to estimate the medium effect alone. It is moreover reassuring 
that the 'Debye' screening length = fi{T = 0) ~ 1.1 fm is also the expected hadronic 
scale. 

At non- vanishing temperature, the screening mass can be determined through a study 
of Polyakov loop correlations {L{0)L~^{r)). Normalization problems make this non-trivial 
1T9| , so that for the moment /i(T) is known only up to an open constant. It is already 



clear, however, that n{T) increases sharply around T = T^, with [djJi/dT) diverging in 
the chiral limit. Combining this result with the known yu(T = 0) and the perturbative 
result /i ~ gT leads to the schematic screening mass form shown in Fig. ^. In QCD, the 
screening mass thus shows a very characteristic behaviour, with a singular derivative at 
T 

The physics leading to this phenomenon seems quite clear. Chiral symmetry restora- 
tion transforms effectively massive into massless quarks. Near T = T^, this sudden loss of 
mass leads to a sudden increase in the density n of coloured constituents and thus, with 
/i ~ n^/^, also in the effectiveness of colour screening. 

An interesting side-line here is the onset of charmonium suppression as signal for 



colour deconfinement Deconfined media will dissociate charmonium and bottonium 
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Figure 8: Schematic temperature dependence of the screening mass /i in full QCD with 
massless quarks of two flavours. 



states in a step- wise pattern, with the larger and more loosely bound states melting before 
smaller and more tightly bound quarkonia. There are indications that the dissociation 



point for the charmonium states Xc(l-P) and ^jj'{2S) coincides with If this could 

be substantiated in more precise lattice studies, it would identify the onset of Xc and 
ip' suppression with the onset of colour deconfinement. 

4. Cluster Percolation 

Conceptually, the deconfinement transition seems rather transparent, no matter what 
the quark mass is. Once the density of constituents becomes so high that several hadrons 
have considerable overlap, there is no mechanism to partition the quark constituents into 
colour-neutral bound states. Instead, there appear clusters much larger than hadrons, 
within which colour is not confined. This suggests that deconfinement is related to cluster 
formation, and since that is the central topic of percolation theory, possible connections 
between percolation and deconfinement were discussed already quite some time ago on a 



rather qualitative level [p2l In the meantime, however, the interrelation of geometric 
cluster percolation and critical behaviour of thermal systems has become much better 
understood and this understanding can be used to clarify the nature of deconfinement. 



To recall the fundamentals of percolation, consider a two-dimensional square lattice 
of linear size L; we randomly place identical objects on of the lattice sites. With 
increasing N, adjacent occupied sites will begin to form growing clusters or islands in the 
sea of empty sites. Define Hp to be the lowest value of the density n = N/ for which on 
the average the origin belongs to a cluster reaching the edge of the lattice. In the limit 
L oo, we then have 

PH~(1--^) , n>np, (28) 



n 



where the percolation strength P{n) denotes the probability that the origin belongs to an 
infinite cluster. Since P{n) = for all n < Up and non-zero for all n > Hp, it constitutes 
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an order parameter for percolation: /3p = 5/36 is the critical exponent which governs the 
vanishing of P{n) at n = Up in two dimensions; in three dimensions, it becomes f3 = 0.41 
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Another quantity of particular interest is the average cluster size S{n), defined as the 
average number of connected occupied sites containing the origin of the lattice; above rtp, 
percolating clusters are excluded in the averaging. This cluster size corresponds to the 
susceptibility in thermal systems and diverges at the percolation point as 

dcc-H ~ |n-npr>, (29) 

with 7p = 43/18 (1.80) as the d = 2 (3) critical exponent for the divergence p5| . 

We now turn once more to the Ising model. For B = 0, the Hamiltonian (0)) has a 
global Z2 invariance (sj —Si V i), and the magnetization m = (s) probes whether this 
invariance is spontaneously broken, as discussed in section 2. Such spontaneous symmetry 
breaking occurs below the Curie point Tc, with 

m(T, 5 = 0) ~ (T, - T)^™ (30) 

governing the vanishing of m(T, B = 0) as T ^ Tc from below. The well-known Onsager 
calculations with /3m = 0.125 for d = 2, a value 10 % smaller than the /3p = 5/36 ~ 0.139 
found for the percolation exponent. 

Since the Ising model also produces clusters on the lattice, consisting of connected 
regions of aligned up or down spins, the relation between its thermal critical behaviour at 
Tc and the onset of geometric percolation is an obvious question which has been studied 
extensively in recent years. In other words, can one interpret magnetization as spin 
domain fusion? This question is now answered 

The geometric clusters in a percolation study consist of connected regions of spins 
pointing in the same direction. In the Ising model, there is a thermal correlation between 
spins on different sites; this vanishes for T 00. Correlated regions in the Ising model 
(we follow the usual notation and call them 'droplets', to distinguish them from geometric 
clusters) thus disappear in the high temperature limit. In contrast, the geometric clusters 
never vanish, since the probability for a finite number of adjacent aligned spins always 
remains finite; it increases with dimension because the number of neighbours does. Hence 
from the point of view of percolation, there are more and bigger clusters than there are 
Ising droplets. 

If percolation is to provide the given thermal critical behaviour, the definition of cluster 
has to be changed such that the modified percolation clusters coincide with the correlated 
Ising droplets |2^, |2^. This is achieved by assigning to pairs of adjacent aligned spins in 



a geometric cluster an additional bond correlation, present with the bonding probability 

Pb = l- exp{-2J/kT}, (31) 

where 2 J corresponds to the energy required for fiipping an aligned into a non-aligned 
spin. The modified 'F-K' percolation clusters now consist of aligned spins which are also 
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bond-connected. Only for T = are all aligned spins bonded; for T > 0, some aligned 
spins in a purely geometric cluster are not bonded and hence do not belong to the modified 
cluster or droplet. This effectively reduces the size of a given geometric cluster or even 
cuts it into several modified clusters. For T — oo, ^ 0, so that the geometric clusters 
still in existence there are not counted as droplets, solving the problem mentioned in the 
previous paragraph. 

For such combined F-K site/bond clusters, full agreement between percolation and 
thermal critical behaviour of the Ising model is achieved for any space dimension d. The 
percolation threshold is now at Tc, the cluster size coincides with that of the correlated 
regions in the Ising model, and numerical simulations show that the critical exponents for 
the new cluster percolation scheme become those of the Ising model. 

Since, as noted, the deconfinement transition in SU{N) gauge theory falls into the 
universality class of the Ising model, it seems natural to look for a percolation formulation 
of deconfinement Jl^, and first studies indicate that this is indeed possible |3^. In 



SU{2) lattice gauge theory, the Polyakov loop constitutes essentially a generalized spin 
variable, pointing either up or down at each spatial lattice site, but with continuously 
varying magnitude. In two space dimensions, this leads to a 'landscape' of hills and 
lakes of various heights and depths. The crucial question in the extension of percolation 
to such a case is how to generalize the bond weight Eq. (|3T|). For a specific lattice 
regularization, the strong coupling limit, it was shown that the action in SU{2) gauge 
theory can effectively be written in terms of nearest neighbour Polyakov loop interactions, 
with {K/4)'^LiLj in place of the Ising form {J/kT)siSj, where k = A/g^ and g denotes the 
coupling in the gauge theory action p4|. We therefore take 



p,^,^. = l-exp{-2(/s:/4)%L,} (32) 

as bond weight between two adjacent Polyakov loops of the same sign. It is known through 
analytic as well as numerical studies that such a form gives the correct bond weight for 
continuous spin Ising models [^, |36[. To test it here, we have carried out studies on a 



number of different lattices for both two and three space dimensions; some results are 
shown in Fig. ^ and in Table 1. Fig. |^ shows that the rescaled percolation probability, 
using the Ising value for the exponent z/ leads to a universal curve, as required. In the 
table we summarize the excellent agreement between thermal and percolation values for 
the critical exponents. The exponents for random site percolation, on the other hand, are 
seen to be considerably different. 

For the specific lattice regularization used, we thus indeed find that deconfinement in 
SU{2) gauge theory can be described as Polyakov loop percolation. Hopefully this can 
be extended to more general lattice regularizations. 

The greatest interest in a percolation approach to deconfinement is, however, based on 
the possibility to define the transition for arbitrary values of the quark mass. Consider the 
case of full QCD with colour S'[/(3) and two massless quark flavours. For niq oo, this 
leads to the first order transition of SU{?>) gauge theory, as counterpart of such a transition 
in the three-state Potts model. For decreasing quark mass, the transition will eventually 
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L-Percolation 
Symmetry Breaking 



Ising Model [33 



Random Percolation 37 



0.528(15) 
0.523(12) 
0.518(7) 
0.4770(10) 



1.985(13) 
1.953(18) 
1.970(11) 
2.0460(39) 



0.632(11) 
0.630(14) 
0.6289(8) 
0.8765(16) 



Table 1: Comparison of percolation and thermal exponents for 3+1 SU{2). 
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Figure 9: Rescaled percolation probability for 3+1 SU{3) gauge theory, using the Ising 
exponent u = 1. 

disappear at a second-order end point defined in terms of m^, and for > > 0, 
there presumably is no thermal transition at all. Finally, for rriq 0, there is the chiral 
symmetry restoration transition. In Fig. |1^, this behaviour is illustrated. Does this mean 
that in the quark mass region > rriq > (which includes our real physical world of 
small but finite bare quark masses) there is no way to define deconfinement as a critical 
phenomenon? 

To address this question, we return to percolation in the Ising model with a non- 
vanishing external field. For H ^ 0, the Ising partition function does not contain any 
singularity as function of T and hence does not show any critical behaviour the Z2 
symmetry responsible for the onset of spontaneous magnetization is now always broken 
and m(T, H ^ 0) ^ for all T. On the other hand, the average size of site/bond clusters 
in the above sense increases with decreasing temperature, and above some critical tem- 
perature it diverges. Hence percolation will occur for any value of 5. In other words, the 
critical behaviour due to percolation persists, while that related to spontaneous symme- 
try breaking and magnetization disappears. At i? = 00, all spins are aligned, leaving the 
bonds as the relevant variables; the system now percolates at the critical density for pure 



15 



oo 







m. 



1 St / order 



confined ! deconfined 



cross 



over 







oo 



Figure 10: QCD phase structure as function of temperature T and quark mass m„ 



bond percolation, which leads to a critical temperature Tk somewhat above the Curie 
point Tc. For finite B, the corresponding values of the critical temperature lie between Tc 
at B = and Tk aX B = oo; they define the so-called Kertesz line p^ , pSj ; see Fig. 0. A 
fundamental and quite general question in statistical physics is what happens at this line. 
Can one generalize critical behaviour to situations where the partition function Z{T) is 
analytic, but where percolation as defined in terms of the input dynamics persists? It is 



evident from the similarity of Figs. and |TT] that the answer is immediately relevant to 
the study of phase transitions in QCD, with deconfinement as the QCD counterpart of 
the Kertesz line. 
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Figure 11: Percolation pattern for the Ising model with external magnetic field B. 



In statistical physics, such generalized critical behaviour indeed leads to observable 
effects. Familiar instances are found in solution/gel transitions, as encountered in the 
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boiling of an egg or the making of pudding. While these phenomena do not result in any 
thermodynamic singularities, they are well-defined transitions which can be quantitatively 
described in terms of percolation [|38[ . 



In closing, we note another possible illustration of the relation between percolation and 
thermal critical behaviour. Consider the F-K site/bond percolation in the two-dimensional 
Ising model as introduced above, and imagine that current can flow between two or more 
bonded sites. In this case, conductivity sets in at the percolation point, independent of 
the magnetization critical behaviour of the Ising model; the system is non-conducting 
below the percolation point and conducting above it. In other words, we now have two 
independent critical phenomena, the onset of conductivity and the critical behaviour of 
Ising thermodynamics, and the former can survive even when the latter is no longer 
present. 

5. Summary 

We have seen that hadronic matter at sufficiently high temperature and low baryon 
density becomes a quark-gluon plasma. In this deconfinement transition, the colour- 
neutral bound-state constituents are dissolved into their coloured components. At high 
baryon density and low temperature, the deconfined medium could be a condensate of 
coloured bosonic diquarks. 

In pure SU{N) gauge theory, colour deconfinement arises through the spontaneous 
breaking of a global Z/^r symmetry of the Lagrangian. In the chiral limit of full QCD, it 
occurs through a strong explicit breaking of this symmetry, due to an effective external 
field setting in when the chiral condensate vanishes. Hence the deconfinement and chiral 
symmetry restoration transitions coincide. 

Colour charge screening in QCD leads to a specific singular behaviour of the screening 
mass. At the critical temperature of chiral symmetry restoration, the effective quark 
mass shift leads to sudden increase in the density of constituents and hence to more 
effective screening. A particularly enticing question here is whether the dissociation of 
the Xc state occurs at just this point - it would then be a measurable order parameter for 
deconfinement . 

Cluster percolation provides an approach to study the geometry of deconfinement. In 
SU{N) gauge theory, first finite temperature lattice calculations indicate that Polyakov 
loop percolation is indeed an equivalent way to identify deconfinement. In full QCD as in 
spin systems with non-vanishing external field, percolation persists even in the absence 
of thermal transitions. It thus seems conceivable to identify the colour deconfinement 
transition for arbitrary quark mass through the onset of percolation. 
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